DUALLY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 03043 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 

THE EFFECTIVENESS OF HEAT EXCHANGERS 
WITH ONE SHELL PASS AND 
FIVE TUBE PASSES 

by 

Lawrence E. Hess 
September 1985 

Thesis Advisor: Allan D. Kraus 

Approved for public release; distribution is unlimited 



T222821 



SECURITY CLASSIFICATION of THIS PAGE (When Dmtm Enter'd) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. REPORT NUMBER 


2. GOVT ACCESSION NO. 


3. RECIPIENT'S CATALOG NUMBER 


4. TITLE (and Subtitle) 

The Effectiveness of Heat Exchangers 
with One Shell Pass and Five Tube 


S. TYPE OF REPORT 8 PERIOD COVERED 

Master's Thesis; 
September 1985 


Passes 




6. PERFORMING ORG. REPORT NUMBER 


7. AUTHORf*; 

Lawrence E. Hess 

V 


8. CONTRACT OR GRANT NUMBERfa; 


9- PERFORMING ORGANIZATION NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 93943-5100 


10. PROGRAM ELEMENT, PROJECT, TASK 
AREA & WORK UNIT NUMBERS 

1 


11. CONTROLLING OFFICE NAME AND ADDRESS 

Naval Postgraduate School 




12. REPORT DATE 

September 1985 


Monterey, California 93943-5100 


13. .NUMBER OF PAGES 


14. MONITORING AGENCY NAME & ADDRESS (if different from Controlling Office) 


IS. SECURITY CLASS, (of this report) 

\ 






Unclassified 

j 






1S« DECLASSIFICATION DOWNGRADING 

SCHEDULE i 


16. DISTRIBUTION STATEMENT (of this Report) 




— 3 


Approved for public release; 


; distribution is unlimited 


17. DISTRIBUTION ST ATEMENT (of the abstract entered in Block 20, It different from Report) 


18. SUPPLEMENTARY NOTES j 

j 


19. KEY WORDS (Continue on reverse aide If neceasmry and Identity by block number) 


I 


Heat Exchanger 
Effectiveness N Method 

Five Tube Pass 






20. ABSTRACT (Continue on reverse eide If necesmmry and Identify by block number) 

Heat exchangers with one shell pass and n tube passes are 
often referred to as 1-n exchangers. The heat transfer 
literature contains many references to studies of 1-n 
exchangers when n is even but apparently, other than a single 
study and some work pertaining to the 1-3 exchanger at the 
Naval Postgraduate School, little work has been done with 
respect to the 1-n exchanger when n is odd. This thesis expands 
upon the study of 1-n exchangers with n being odd. 



DD I j an ^73 1473 EDITION OF I NOV 65 IS OBSOLETE 



S N 0102- LF- 014- 6601 



1 SECURITY CLASSIFICATION OF THIS PAGE Dmtm Entmrmd) 



SECURITY CLASSIFICATION OF THIS PAGE fW h«n Data Enffd) 



20. (Continued) 

While a completely closed form solution was found to be 
unfeasible, a polynomial approximation has been developed 
that yields the effectiveness (e) of the two possible arrange- 
ments of the 1-5 exchanger as a function of the capacity rate 
ratio (R) and the number of transfer units (N t ) • This will 
enable the analyst to consider exchangers where the inlet 
to and outlet from the tubes are at opposite ends of the 
exchanger . 



5 N 0102- LF- 014- 6601 



2 



SECURITY CLASSIFICATION OF THIS PAGEfRTian Dal* Enffd) 



Approved for public release; distribution is unlimited 



The Effectiveness of Heat Exchangers 
With One Shell Pass and 



Five 


Tube Passes 




by 



Lawrence E. Hess 

Lieutenant Commander, United States Navy 
B.S., The Citadel, 1973 

Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE 


IN MECHANICAL ENGINEERING 




from the 



NAVAL POSTGRADUATE SCHOOL 
September 1985 



ABSTRACT 



7 * - ; i 
- 2 - 

f*»/ 



Heat exchangers with one shell pass and n tube passes 
are often referred to as 1-n exchangers. The heat transfer 
literature contains many references to studies of 1-n 
exchangers when n is even but apparently, other than a 
single study and some work pertaining to the 1-3 exchanger 
at the Naval Postgraduate School, little work has been done 
with respect to the 1-n exchanger when n is odd. This 
thesis expands upon the study of 1-n exchangers with n being 
odd . 

While a completely closed form solution was found to be 
unfeasible, a polynomial approximation has been developed 
that yields the effectiveness (e) of the two possible 
arrangements of the 1-5 exchanger as a function of the 
capacity rate ratio (R) and the number of transfer units 
(N tu ). This will enable the analyst to consider exchangers 
where the inlet to and outlet from the tubes are at opposite 
ends of the exchanger. 
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NOMENCLATURE 



English Letter Symbols 

A = Exchanger heat-transfer surface, sq m 

A m = Coefficient of the rath value dimensionless 

Aq = Coefficient, dimensionless 

- 1 st order coefficient to be multiplied by N^ u , 
dimensionless 

A 2 = 2 nd order coefficient to be multiplied by N^ u 2 , 
dimensionless 

A 3 = 3rd order coefficient to be multiplied by N-^3, 
dimensionless 

Aq = 4th order coefficient to be multiplied by N^ u ^, 
dimensionless 

A 5 = 5th order coefficient to be multiplied by N^ u ^, 
dimensionless 

a = Exchanger heat-transfer surface, sq m/m 

C = Capacity rate, W/°K. Also designates dimensionless 
arbitrary constant 

Cp C = Specific heat at constant pressure of cold fluid, 
J/kg°K 

C D v, = Specific heat at constant pressure of hot fluid, 

J/kg°K 

D = Empirical value of effectiveness (computer generated) , 
dimensionless 

e = Error. Also used as the exponential function 

F = Logarithmic mean temperature difference correction 
factor, dimensionless 

L = Exchanger length, m 

N = Number of effectiveness empirical data points used to 
determine a curve for R, dimensionless 
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n = Number of tube passes, dimensionless. Also number of 
equations, dimensionless 

n c = A related N-t u per unit length cold side, m”l 
n^ = A related Nt u per unit length hot side, m“l 
N tu = Number of transfer units, dimensionless 
P = Temperature group, dimensionless 

q = Total rate of heat transfer, W 

Hmax = Maximum total rate of heat transfer, W 
R = Capacity rate ratio, dimensionless 

S = Temperature group, dimensionless 



S r 


= 


Sum 


of the 


squares of 


the residuals, dimensionless 


T 


= 


Hot 


fluid 


temperature , 


°C 




T pi 


= 


Particular 


integral, dimensionless 




Tl 


= 


Hot 


fluid 


temperature 


in, °C 




t 2 


= 


Hot 


fluid 


temperature 


out, °C 




tl 


= 


Cold 


fluid 


temperature 


in, °C 




t2 


= 


Cold 


fluid 


temperature 


out, °C 






= 


Cold 


fluid 


temperature 


1st pass, °C 




tab 


= 


Cold 


fluid 


temperature 


between 1st and 2nd 


passes 


tb 


= 


Cold 


fluid 


temperature 


2nd pass, °C 




tbc 


= 


Cold 


fluid 


temperature 


between 2nd and 3rd 


passes 


tc 


= 


Cold 


fluid 


temperature 


3rd pass, °C 




ted 


= 


Cold 


fluid 


temperature 


between 3rd and 4th 


passes 


td 


= 


Cold 


fluid 


temperature 


4th pass, °C 




tde 


= 


Cold 


fluid 


temperature 


between 4th and 5th 


passes 


te 


= 


Cold 


fluid 


temperature 


5th pass, °C 
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U = Overall heat transfer coefficient, W/m^ - °C 
W = Mass flow, kg/sec. 

x = length coordinate, m. Also used to represent 
a constant value in a sequence 

y = Sum of mth degree polynomial, defined by eq. (5.1), 
dimensionless 

Greek Letter Symbols 

a = Root of auxiliary differential equation, 1/m 

e = Exchanger effectiveness, dimensionless 

X = Combination of variables defined by equation (3.12), 
dimensionless 

1 = Summation, dimensionless 

a2 = Variance, dimensionless 

8 m = M ean temperature difference for exchanger, °C 
3 = Indicates partial derivative, dimensionless 

0 = A combination of terms defined by eq . (3.23) 

Subscripts 

c = Cold fluid 

h = Hot fluid 

i,j,k = Values in a sequence 

m = Degree or order, an exponent 

1 = inlet 

2 = outlet 
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I. INTRODUCTION 



In many engineering applications, it is desirable to 
transfer thermal energy between fluids at different 
temperatures. This action may be economically accomplished 
with the aid of a heat exchanger. 

The primary function of a non-direct contact heat 
exchanger is to provide two paths of flow, one for the hot 
fluid and one for the cold, by means of which heat can be 
transferred through walls separating the fluids. 

The quantity of heat transferred is governed by three 
factors: 1) the amount and nature of the heat transfer 

surface to which the two fluids are exposed, 2) the overall 
heat transfer coefficient, and 3) the mean temperature 
difference across the intervening wall from one fluid to the 
other . 

Numerous flow arrangements and geometric designs are 
possible. One of the most common types of heat exchangers 
is the shell and tube type. Here one fluid flows through 
the tubes. The tubes are enclosed in a shell with 
provisions for the other fluid to flow through the spaces 
between and exterior to the tubes. Overall flow may be 
either parallel or counter, but in most cases, the latter is 
preferred. The reason for this will be seen later. One of 
the most common flow arrangements used in the shell and tube 
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heat exchangers is a combination of the two, the parallel 
counterflow exchanger. Heat exchangers with one shell pass 
and n tube passes are therefore referred to as 1-n heat 
exchangers . 

One example of this flow arrangement is the one 
shell-pass two tube pass exchanger as shown in Figure 1.1. 




Figure 1.1 A 1-2 Heat Exchanger (One Shell Pass and Two 
Tube Passes) 

This type of exchanger is configured so that all of the tube 
side fluid flows consecutively through half of the total 
number of tubes: The baffles in the shell serve to mix 

the shell fluid and to maintain fluid flow at right angles 
to the tube passes in order to obtain a maximum shell 
side" heat transfer coefficient. 
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O'Hare [Ref. 1] has indicated that, to date, much work 
has been done on finding the logarithmic mean temperature 
difference for heat exchangers with an even number of passes 
but little has been done with exchangers having an odd 
number of passes. From this basic statement, he proceeded 
with a theoretical examination of the 1-3 heat exchanger 
which resulted in some startling conclusions. As a follow 
up to his work, this thesis will parallel and follow the 
format of O'Hare [Ref. 1] for the 1-5 parallel counter flow 
exchanger. One arrangement of the 1-5 exchanger is 
displayed in Figure 2. 



l i T 1 ^Baffle 




Figure 1.2 A 1-5 Heat Exchanger (One Shell Pass and Five 
Tube Passes) 
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II. THE DEVELOPMENT OF THE EFFECTIVENESS METHOD 



A. HISTORICAL DEVELOPMENT 

Nagle [Ref. 2: pp. 604-609], in 1931, credited Davis 

[Ref. 3] with a simplified method for computing actual temp- 
erature differences between two heat-interchanging streams 
which depart from true counter or concurrent (parallel) 
flow. This is now the familiar "F factor" method which 
expresses the actual mean temperature difference 6 ra in q = 
UA0 m as a fraction F of the counterflow logarithmic mean 
temperature difference, LTMD , 8 mc via 9 m = F0 mc . 

The example of initial interest was the 1-2 exchanger 
with a single shell pass and two continuous tube passes in 
counter and concurrent flow with it. The method involved 
derivation of the actual temperature difference for the flow 
pattern and formed the ratio F = 9 m /0 mc . This familiar LMTD 
correction factor was plotted conveniently as functions of 
the effectiveness, e, and the capacity rate ratio R with R 
as a parameter. These mean temperature difference 
correction charts are available for many flow arrangements 
[Ref. 4: pp. 829-833 and Ref. 5]. The effectiveness, e 

(often called P or S) , is always the cold fluid 
effectiveness and R is always the capacity rate ratio of 
cold fluid to hot fluid. 
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Nagle detailed assumptions and derivations for the 1-2, 
1-4 and 1-6 exchangers. The F factors were obtained by 
Nagle through graphical integration and were accompanied by 
the comment that F factors for the 1-2 exchanger could be 
applied with negligible error to 1-4 and 1-6 exchangers. 
Underwood [Ref. 6: pp. 145-148] rederived the equations of 
Nagle for 1-2 and 1-4 exchangers to eliminate the need for 
obtaining F factors by graphical integration. 

Bowman [Ref. 7: pp. 541-544] pointed out that for a 
very large or infinite number of tube passes, the F factor 
approached, as a limit, its value in crossflow with both 
fluids completely mixed. It was further stated that even at 
the limit, the F factors were only 1 to 2 percent lower than 
those for the 1-2 exchanger. A previous paper by Kraus and 
Kern [Ref. 8] did not confirm the generalization that 1-n 
exchangers differed only negligibly from the 1-2 exchanger 
although this lack of confirmation was obtained on an 
e = f(R, N^- u ) basis. Moreover, the Kraus-Kern work does not 
confirm the generalizations on an F = f(R, N-£ U , e) basis. 

From the standpoint of usefulness and good accuracy, it 
is essential that F factors, if they are to be used in 
preference to e = f(R, N-^u , flow arrangement), be obtained 
with precision. Plots of F = f(R, e, = P or S) [Ref. 4: 
pp. 829-833 and Ref. 5] show that the curves for particular 
values of R approach infinite slope as F decreases. While 
this can be partially alleviated by restricting R < 1.0 (a 
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constraint used in the e = f(R, N tu , flow arrangement) 
approach), it is seen that small errors in the interpolation 
for R or e = P or S can result in large fluctuations in the 
value of F. 

In a comprehensive paper, Bowman, Mueller and Nagle 
[Ref. 9: pp. 283-294] presented graphs of F factors for 
shells with one through six shell passes and numbers of 
continuous tube passes respectively double the number of 
shell passes. In view of the earlier references to Nagle 
and Bowman, it should be noted that F factors were computed 
for the 1-2 exchanger in [Ref. 9: pp. 283-294] using the 
equations of Underwood [Ref. 6: pp. 145-158]. 

Ten Broeck [Ref. 10: pp. 1041-1042] prepared a graph 
of the dimensionless groups now known as e , R and N tu for 
the 1-2 exchanger. Such a graph had the added versatility 
of simplifying the calculation of performance in a given 
exchanger when operating at conditions different for those 
for which it was designed. Kays and London [Ref. 11: 
pp. 63-74] prepared similar graphs and tables of 
s = f(R, Nj- U , flow arrangement) for the 1-2 exchanger and 
for several cases of crossflow and periodic flow. 

In the foregoing discussion, it has been observed that 
there are two methods for the design and analysis of a heat 
exchanger. These are the so-called "F factor" and "e-N^- u " 
approaches and it is to be noted that the Heat Transfer 
Literature contains many references of design equations for 
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either of these methods for exchangers that have one shell 
pass and even numbers of tube passes. 

It was O' Hare [Ref. 1] who pointed out that there may be 
a virtue in having an odd number of tube passes particularly 
in a marine application. Such an arrangement allows the 
tube side fluid to enter one side of the exchanger and leave 
from the opposite side (see Figure 1.2). 

With the possible use for an odd tube pass heat 
exchanger established, a literature search was conducted to 
uncover any previously solved cases of 1-n heat exchangers 
with n, the number of odd tube passes. O' Hare [Ref. 1], 
Fischer [Ref. 12], and Stevens, et. al. [Ref. 13] have 
dealt with this type of exchanger. Although Stevens, 
et . al . deal with multiple crossflow passes, all provide the 
foundation for the work reported here as well as additional 
work pertaining to the 1-3 exchanger. O'Hare's investiga- 
tion centered around the development of a solution for the 
1-3 exchanger in the ef f ectiveness-N tu framework while 
Fischer derived an equation using the mean temperature 
difference for the 1-3 exchanger. Moreover, Fischer 
expressed his results in terras of the correction factor F 
and treats only the case of the 1-3 exchanger. This thesis 
will expand on the N tu method of thermal heat exchanger 
design by covering the 1-5 heat exchanger. The work will 
treat both the three counterflow and two parallel flow 
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and the two counterflow and three parallel flow (1-5: 3P) 
cases (see Figures 2.1 and 2.2). 
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Figure 2.2 1-5 :3P Five Tube Passes - Three in Parallel Flow 
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III. AN ATTEMPT AT AN ANALYTIC SOLUTION 



A. EFFECTIVENESS AS A FUNCTION OF CAPACITY RATE RATIO AND 
EXCHANGER SIZE (N tu ) 

This section deals with an investigation into the 
effectiveness, e , of a one shell pass-five tube pass heat 
exchanger, in which e compares the actual heat transfer 
capability to the thermodynamically limited, maximum possible 
heat transfer capability. This exchanger heat transfer 
effectiveness is given by 

_ q _ ^h^ T hot,in T hot,out^ _ ^c ^cold , out ^coldjin^ 

q max min^ not,in cold, in' min^ hot, in cold, in) 

where C rain is the smaller of the C^ and C c magnitudes (the 
capacity rates). Thus, e possesses the significance of 
effectiveness of the heat exchanger from a thermodynamic 
point of view, with the magnitude of the effectiveness 
completely defining the heat transfer performance. In 
general we express e = f(N tu , R, and flow arrangement) and 
when the flow arrangement is specified or understood, it is 
said that e = f(N^ u , R). [Ref. 11: pp. 14-26]. 

The number of heat transfer units N tu is a nondimensional 
expression of the "heat transfer size" of the exchanger. 

When N^ u is small the exchanger effectiveness is low, and 
when N tu is large, e asymptotically approaches the limit 
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imposed by the flow arrangement and thermodynamic 
conditions- From inspection of the definition of N tu 





1 



A 

j UdA 
0 



it is clear that the overall conductance and transfer 
area affect the costs of attaining a high value for N tu , 
ergo a high e. The capacity rate ratio, R, as defined by 



is simply the ratio of mass flow rate times specific heat 
capacity for the two streams. These can be considered as 
flow stream thermal-capacity rates, i-e- , the energy storage 
rate in the stream per unit of temperature change- [Ref. 

11: pp. 14-26] 

The attempt taken in this thesis to develop a closed 
form solution has used the basic fundamentals of heat 
transfer as well as the foregoing definition. A closed form 
solution for e was sought for both 1-5 exchangers with one 
having three out of five tube passes in parallel flow and 
the other having three out of five tube passes in counter- 
flow. The analytical approach taken, and demonstrated in 
this section, is for three out of five tube passes in 
counterflow . 







max 
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B. ANALYTICAL DEVELOPMENT 



The derivation for the effectiveness, e, of the 1-5 
exchanger as a function of the capacity rate ratio, R, and 



number 


of transfer units, N^ u , depends on several 




assumptions . 




(1) 


The overall coefficient of heat transfer, U, 
not vary within the exchanger. 


does 


(2) 


The specific heat of both hot side and cold 
fluids does not vary. 


side 


(3) 


Each fluid is thoroughly mixed, that is, the 
temperature of both hot and. cold side fluids 
uniform over any cross section. 


is 


(4) 


Steady flow conditions are maintained. 




(5) 


Heat losses to or from the environment are 
negligible. 




(6) 


No change of phase takes place; all heat 
transferred is sensible heat. 




(7) 


There is equal heat transfer surface in each 


pass . 



The configuration is shown in Figure 3. 1 where the five 
tube passes are designated with subscripts a, b, c, d and e. 

The temperature of the hot (shell side) fluid is indicated 

by upper case letters. For the cold (tube side) fluid, 
lower case letters are used. The subscript 1 always refers 

to the fluid inlet and the subscript 2 always refers to the 

fluid outlet. 

With Wfo and Cp^ designating mass flow (kg/sec) and 
specific heat (Joules/kg -<> K) of hot fluid entering at T^ and 
leaving at T 2 we define a capacity rate for the hot side 



^h - WhCph 
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Figure 3.1 Three Counter Two Parallel Pass Configuration 
for Development of the Sought After 
Effectiveness Relationship 
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In similar fashion for the cold side (with W c and C pc ) 
entering at t^ and leaving at t 2 , we have 

C c = ^c^pc 

We then obtain an energy balance for the entire 
exchanger 

C h (Ti - T 2 ) = C c (t 2 - ti) (3.1) 

Over the right hand side of the exchanger (Figure 3.1) 
C h (Ti - T) = C c (t 2 - t e + t d - t c + t b - t a ) (3.2) 
and a differentiation gives 

C h dT = C c (dt e " dt d + dt c - dt b + dt a ) (3.3) 

Across dx, with a (m^/m) , the surface per running meter 
of length of pass so that A = 5aL is the total surface in 
the exchanger, we may write the heat transferred to the 
element dx in each cold pass. 



C c dt a = 


Uadx( T 


^a ) 


(3.4a) 




-Uadx(T 


“ t b ) 


(3.4b) 


C c dt c = 


Uadx( T 


- t c ) 


(3.4c) 


= 


-Uadx( T 


- t d ) 


(3 . 4d) 


= 


Uadx( T 


— t e ) 


(3. 4e) 



26 



Here it should be observed that due cognizance has been 
taken of the direction of the flow in each cold fluid pass 
with respect to the positive sense of the length coordinate, 
x, and U is the overall heat transfer coefficient (W/m 2 -°C). 
With eqs. (3.4) in eq. (3.3) 

C^dT = Ua(5T — t^ — t^ — t^ t^-j — t e )dx 

or 

3 I = n h (5T - ( a - l b - ‘o' ' *e) < 3 ' 5 > 

where 



n _ Ua 

"h - 



is a sort of N^ u per unit length for the hot side. 
Now differentiate eq. (3.5) 



d 2 T dT _ ^a _ ^b _ ^c _ ^d _ ^e . 

, 2 n h^dx dx dx dx dx dx 

dx 



and with eqs. (3.4) substituted 



d 2 T 



dx" 



= 5n 



dT 



h n c n h (T t a + t b t c + t d t e )(3,6) 



where 




where again the resemblence of n c to N^ u can be noted. 
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From eq. (3.2) we obtain 



C~~ < T 1 - T) - t 2 = t b - t a - t c + t d - t e 
c 



(3.7) 



and with eq. (3.7) put into eq. (3.6) 



d 2 T , dT 

— 2 - 5n h Sk 
dx 



- “c“h l T + JT < T 1 - T > - V 
c 



or 



d 2 T 



dT 



dx 



“ 5n h Hi ' " n c n h TT [(R c 1)T + T 1 R c t 2 1 

c 



where 

R c = C c /Ch 

is the capacity rate ratio for the cold side. 
Notice that 



C h Ua Ua C h ,Ua 2 

n c n h C~ ‘ ' C~ ' ~ ’ ( C-> 01 

c c h c c 



and 



i C. 

r - 1 - h 

n h R C 

c c 



(3.8) 



a capacity rate ratio for the hot side. Then, algebraic 
adjustment provides 
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(3.9) 



d 2 T _ dT , 

— 2 ' 5n h cE + m < 
dx 



1 - R 



RT Ji > T - “<1C - V 



which is a linear, non-homogeneous , second order 
differential equation with constant coefficients having a 
complementary function 



a. x a 9 x 

T c - C l e + C 2 e 



(3.10) 



where and C 2 are arbitrary constants and where 



5n, . 9 1 - R. 1/2 

a i ’ a 2 ~r~ ± 2 [25n h " 4ra(— i;- — )] 



h 



But 



5n. n, A 1 - R. 1/2 

h ± h [25 . 4m ( _h , 



2 " 2 



2 v R. 
n h h 



/TT .2 (C, ) C, 2 _ . 

m _ (Ua) # v h J _ , h s _ ™ 2 _ 1 



,"2 ’ 77772 " ( C _) ' R h 

n h ( C c ) ( Ua > c 



R 



so that 



a l ’ a 2 



5n. n. 

h t — [25 



T" 



M 2 , 1 - V // 2 

4R h <— K—>l 



or 



“l ,a 2 



(5 ± X) 



(3.11) 
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where 



A = [25 - 4R h (l - R h )] 1/2 (3.12) 

Designate the particular integral as T P1 and by the 
method of undetermined coefficients let T p ^ = P so that in 
eq. (9) 



1 - 



»h 



m ( — x ) P = m(j5 T n ) 



R h ^ 



This makes 



T . = P 
Pi 



to R h 

^ R“ ~ T l^ U - bJ 
h h 



so that 



*2 - R h T l 



pi = 



1 - 






(3. 13) 



The general solution to eq . (3.9) is the sum of eqs. (3.10) 

and (3.13) 

ct . x Ct 9 x t 9 - R. T 1 

T(x) = C ie + C 2 e Z + ~ r - ="g- - (3-14) 

h 



where the arbitary constants, and C 2 are evaluated from 
conditions at x = 0 and x = L. At x = 0, T(x =0) = T 2 and 
at x = L, T(x = L) = Tj_. When these are inserted, in turn, 
into eq. (3.14), one obtains a pair of linear algebraic 
equations in the unknowns and C 2 
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T 2 ’ C 1 + °2 * T pi 



a-L a^L 

T. = C.e + C 0 e + T . 

11 2 pi 



where Tp^ is given by eq. (3.13). 

It is only a matter of algebra to show that 



C 



1 



Ctr>L 

(T- - T . ) - (T 0 - T . )e 

v 1 pi _ _ 2 pi 1 

a^L oi^C 

e - e 



and 



C„ = 



a i L 

(T_ - T . )e - (T. - T . . 

v 2 pi 7 _1 pi) 

a^L a 2^ 
e - e 



It is easy to see from eq. (3.1) that 



D C h (t 2 t l ) 

h ^ ( T i - Tp 



so that 



t 2 - ti + R h (Ti - T 2 ) 



Use of this in eq. (3.13) shows that 



T . = 
Pi 



h + V T i " V - R h T l 



1 - 



*h 



(3.15a) 



(3.15b) 
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or 



h - v 2 

Pi - 1 - H„ 



(3.16) 



indicating two alternative forms for Tpi given by eqs. 
(3.13) and (3.16). 

Insertion of eqs. (3.13) and (3.16) in eqs. (3.15) for 
Ci and C 2 will yield after some algebra 



C 



1 






T 2 a 2 L 

a 2 h 

e 



(3.17a) 



and 



C 



2 




t l a l L T 
— ) e - (- 

h 



a^L a 2^* 
e - e 




(3.17b) 



Equation (3.14) is an expression for the hot side 
temperature at any location in the exchanger in terms of the 
extreme temperatures, t^ , t 2 , Tj_ and T 2 . 

Next take eq. (3.5) and set it equal to the derivative 
of eq. (3.14) noting that C]_, C 2 and T p ^ are all known 
constants . 



CX X Ot X 

dx = n h^ 5T “ t a “ t b “ t c “ td ” te ^ = a i c l e + a 2 C 2 e 

(3. IS) 

At x = 0, where T = T 2 , t a = ti, tfc = t c = tb c , = t e = t^e 
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(3.19) 



dT 



dx n h^ 5T 2 t l 2t bc“ 2tde ^ a l C l + a 2 C : 



and if we subtract eq. (3.4a) from eq. (3.4c) we obtain 



dt - dt 
a c 

t - t 
a c 



Ua , , 

7T- dx = -n dx 

v ^ C 



which can be integrated using C3 as the constant of 
integration 



- ‘c = C 3 e 



-n x 
c 



At x = 0 where t a = ti and t c = t^ c 



^1 “ ^bc ^3 



or 



t bc = tl “ c 3 



In addition at x = L, t a = t a b and t c - t C( j so that 



t , - t , = C„e 

ab cd 3 



-n L 

c 



or 



t . - t , n L 

r - ab cd = ( _ t )0 c 

C 3 =n~L C ab cd je 



This gives a relationship between t a b and t^ c 



N 



t bc “ t l (t ab " t cd )e 



(3 



. 20 ) 
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where N c = n c L can be considered as the total number of 
transfer units for the cold side. 

The same procedure may be employed to find t^g. 
Subtract eq. (3.4c) from eq. (3.4e) 



dt 

~t 



c 



c 




Ua 

XT 

c 



dx 



n dx 
c 



and by integration 



t 



c 



t 

e 



-n 



C 



4 



e 



c 



x 



At x = 0 where t c = t^g and t e = t^e 
^bc “ ^de = C4 



or 



tde - ^bc “ ^4 



At x = L where t c = t cc j and t e = t2 



-n c L 

^cd “ ^2 ^4 e 



and again with n c L = N c 



N c 

C4 — (t cc } t2)e 

and 

N c 

^de “ ^bc (^cd ^2) e 



(3.21) 
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As before substract eq. (3.4b) from eq. (3.4d) 



dt 

~t 



b 

b 





n dx 
c 



and by integrating find that 



*b - 



C 5 e 



n c x 



at x = 0, t b = t bc and t d = t c j e so that 
^bc " ^de = ^5 

and at x Li y t b t ab ^^d t d t^> d * This shows that 

n c L 

tab “ t cd = ^ 5 e 

and 

N c 

tcd “ t ab (t b c t de )e (3.22) 

where N c = n c L. 

Return to eq. (3.18) and look at the conditions at x = L 
where t a = t b = t ab , t c = t d = t cd , t e = t 2 and T = Tj.. 

These conditions in eq. (3.18) give 



V 5T 1 - 



2t 



ab 



- 2t 



cd 



t 2 ) a i c i e 



a 



a 2 C 2 e 



cipL 



(3.23) 

where again we remember that and C 2 are known constants. 

Now with four equations, eq. (3.20), eq. (3.21), 
eq. (3.22), and eq. (3.23) and four unknown variables we 
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can solve for the four variables, t a b> tb c > t C( j, and t e . 
Cramer's rule will be used to solve for two of them; t^c and 
tde • 



First arrange eqs. (3.20) to (3.23) in standard form 
with the unknowns appearing in columns. 



~ N c 

^b + tbc e " 


t , 
cd 


-N c 

= t t e 


(3.20) 


’ Sc + 


N c 

t ,e + t , 

cd de 


N c 

= V 


(3.21) 


tab “ ^bc + 


ted + tde 


= 0 


(3.22) 


tab + 


ted 


= 0 


(3.23) 



where 0 



1/2(5T 1 




( a ^ ® 



L a 2 L 

a 2*"2 e ^ ) 



If 



AX = B = 



1 


0 

1 

CD 


-l 


0 




I i 

d 

1 




-N c 

t x e 


0 


-l 


e N = 


1 




tbc 




N c 

t 2 e 


1 


O 

CD 


-1 


e N = 




t cd 




0 


1 


0 


1 


0 




Jldg^ 




3 



The solution for tb c and t^e is 



det A. 



be 



'be 



det A 



det A 



de 



'de 



det A 
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where, A bc and A^g are the matrices obtained by replacing 
the entries in the second and fourth columns of A by the 
entries in matrix B. 

Compute the determinant, det A of matrix A (the matrix 
of coefficients) 



det A 



1 



e 



-N 



c 



1 



0 




10 10 



det A 





) 



Since det A * 0, the system should have a unique solu- 
tion. Continue by computing det A bc and det A b e . Thus 



det A bc = 



~ N c 

1 t x e c -1 0 



N r N c 

0 t 2 e c e c 1 



0 -1 



N, 



2N„ N r -Np 

= e c (3 - 2t 9 ) - t. (e c + 2e c ) 
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and 



-N, 



-1 t^e 



-N r 



det A de = 



0 



-1 

Nr 



1 -e 



Nr 



-1 



t„e 



0 



Nr 



0 



2N, 



N C , o " N c. 



= (e ^ + 1) ( g - 2t 9 ) - t. (e ^ + 2e 



Therefore 



t 



be 



det A, 
be 

det A 



2Np N r -Np 

e c (g - 2t 2 ) - t 1 (e c + 2e c ) 



-(e c + 2e c ) 



or 



and 

2Np N c -Np 

det A de (e c + l)(g - 2t 2 ) - t^e c + 2e c ) 

t de = det A = N =N 

-(e c + 2e c ) 



Nr 



'be 



t l 



(3 - 2t 9 ) 



-2N 



(3.24) 



1 + 2e 



or 



t 



de 



2Np 

(e c + 1 ) ( 3 

N =N 

e c + 2e 



c 



2t 2 ) 
— ) 



(3.25) 
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Then with eq. (3.24) and (3.25) in eq. (3.19) and some 
algebraic manipulation 

N c -N r 

(3 - 2t„)(2e + e c ) 

a l C l + a 2 C 2 = n h^ T 2 " t l^ ) + =2N ^ 

(1 + 2e C ) 

(3.26) 

Equation (3.26) is similar to the one derived by O'Hare 
[Ref. 1: eq. (22)]. At this point he continued to attempt a 
closed form solution for e as a function of R and N tu . It 
is to be noted, however, that he was unable to obtain the 
solution because of the presence of terras deriving from 
dimensional parameters such as n^ with dimensions m“l 
[Ref. 1: pp. 46]. Since n^ is present in eq. (3.26), it 
must be concluded that a further attempt at an analytic 
solution for the effectiveness of the 1-5 exchanger will 
suffer a similar fate and be also unachieveable. 

This conclusion should be obvious because, when dealing 
with an equation derived from n equations with n+1 unknowns, 
one cannot create an additional equation by merely 
multiplying one of the n equations by a constant. This 
leads to a linearly dependent set which has no unique 
solution . 

Thus, after observing that O'Hare's attempt failed, this 
thesis, like O'Hare's, turns to a numerical analysis to 
achieve the requisite objective. 
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IV. NUMERICAL AND COMPUTER ANALYSIS 



As O' Hare discovered with the 1-3 heat exchanger 
[Ref. 1], an analytic solution for effectiveness as a 
function of the capacity rate ratio, R, and N tu was also 
unattainable for the 1-5 exchanger. Hence, an alternate 
approach had to be taken and this involved utilizing a 
finite-difference thermal analyzer combined with a 
linearizing scheme for the solution of a set of nonlinear 
algebraic equations [Ref. 14]. This thermal analyzer has 
been employed to provide solutions for temperature that lead 
to effectiveness, e, value for the 1-5 :3C and 1-5: 3P heat 
exchangers as a function of R and N^ u . 

A. THERMAL ANALYZER TVSSI 

The computer program TVSSI used in the development of 
the e = f(R, N^ u ) graphs for the 1-5 heat exchanger [Ref. 1: 
Appendix A], was adapted from TVSS2 listed by Kern and Kraus 
[Ref. 14: Appendix C]. This adaptation consisted of modi- 

fying TVSS2 to compute solutions in the SI system and accept 
an input file specifically designed for the 1-5 heat exchan- 
ger. This program also had to be adapted to the IBM 3033 AP 
system. The program utilizes the Cholesky decomposition 
scheme, and because of the linearization of the radiation 
terras (not used in this study) the program is iterative. 
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B. DEVELOPMENT OF THE 1-5 :3C and 1-5 :3P HEAT-EXCHANGER 

MODELS 

As indicated in the foregoing, TVSSI uses an input file 
to generate a solution for temperatures that will lead to 
the sought after effectiveness as a function of R and N^ u . 
This input file is created by a program which requires an 
input of given capacity rates (hot and cold), overall 
coefficient, surface area and inlet stream temperatures (hot 
and cold) values. 

This required the development of a program which created 
such a file. O'Hare [Ref. 1: pp. 53-55] first developed a 
program using the 1-4 heat exchanger that would create an 
input file for TVSSI. He compared the effectiveness from 
the computer generated results and a known analytic solution 
of effectiveness developed by Kraus and Kern [Ref. 8]. Once 
O'Hare established confidence in his results, he proceeded 
to develop a program which would create an input file for 
the 1-3 exchanger. In this study, a similar technique was 
used to develop a computer program for the 1-5: 3C and 1-5: 3P 
heat exchangers. These programs are called NTU53C and 
NTU53P (see Appendices A and B) and, in them, the following 
parameters and techniques were utilized. 

1. 300 nodes were used. 

2. The initial temperature for the computer to begin the 
iterative process was set at 200°C. 

3. An eventual accuracy of .05 between the final and next 
to last iterations was used. 
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4. A radiation coefficient convergence factor of 0.66667 
between iterations was used. This is required by the 
program even if radiation is not considered as a 
node . 

5. The maximum number of iterations that the computer 
was allowed to perform was set at 12. 

6. A damping factor of .8 was set as an initial damping 
based on the number of non-linear terras in all of the 
node equations. 

With the values of C^, C c , U, A, T^ and ti set, the 
input is set and an input file for TVSSI was generated. In 
this file all node equations and internode conductance 
values were determined. -This program specified the nodes 
that interact with each other and the methods by which the 
thermal interaction takes place such as in the forced 
convection and fluid flow modes. 

Both programs make use of the fact that each term 
in a node equation shows three things. The first is the 
node that is coupled for heat flow with the node in 
question. The second is the method of heat flow between the 
nodes. In this case, forced convection and fluid flow are 
used. Finally, the node equation shows the magnitude of the 
internode heat flow. Here all the pieces of information are 
collected and presented for use by TVSSI as an input file 
with all items in the proper format. 

It is the output values of both NTU53C and NTU53P (a 
typical result is shown in Appendix C) that are used by 
the thermal analyzer (as an input file) to determine 
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the temperatures T 2 and t2 for the specific set of given 
initial parameters C^, C c , R, A, Tj., and t]_. 

C. SCOPE OF COMPUTER ANALYSIS 

With the development of NTU53C and NTU53P computer 
programs, the files that they generated could then be used 
to solve for temperatures which in turn could yield effec- 
tiveness values for a given R and N tu . Of course, it must 
be realized that many computer runs are required to generate 
enough data to ensure confidence in the results which cover 
a wide range of capacity rate ratios and N tu values. 

To efficiently expedite the computer task, the Multiple 
Virtual System (MVS) with Job Entry Subsystem and Networking 
(JES3) was employed. The MVS coupled with JES3 is more 
commonly referred to as batch processing. Based on trial 
and error, it was determined that in order to build a solid 
data base, eleven different values for effectiveness were 
needed to best represent a particular value of R. Thus, a 
range of R from R = 0.01 to 1.0 was used. Due to the 
enormous amount of computer time required to solve a 300 
node problem, R increments of 0.05 were used for this study 
and this thesis reports on this basis. In all, 42 curves for 
both the 1-5: 3C and 1-5: 3P heat exchangers were required. 
With this change in the increments between R values 
(O' Hare's work showed that an increment in R of 0.01 was not 
necessary), a linear interpolation was found to be accepta- 
ble to less than 0.02 percent error. 
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TVSSI was slightly modified in accordance with the 
appropriate guidelines of the job control language (JCL) 
needed to run on the batch processing system. These 
modifications were few and were needed only at the beginning 
and end of TVSSI. The modified version of TVSSI has been 
called TVCOUNT with changes shown in [Ref. 1: Appendix F]. 

It was TVCOUNT that was then used to activate TVSSI. 

It also became necessary to modify the two input file 
programs NTU53P and NTU53C so that they needed to be 
compiled only once. They were then loaded in a library 
file to be used when called by another program. New 
programs utilizing the batch system were written that could 
easily be loaded with the appropriate input data for a 
specific R value. These, which are referred to as "sister 
programs," were used to go from the library file to TVSSI 
and cause TVSSI to be executed eleven times under TVCOUNT 
covering the desired range of N tu for a specific R value. 

The modified input files and their associated "sister 
execution programs are found in Appendices D through G. 

A flow chart of how all of the foregoing is accomplished 
is provided by O'Hare [Ref. 1: Figures 4.4 and 4.5]. It is 

noted that in these figures, TVSSI is referred to as TVSSIA 
through TVSSIV. These are the same programs as TVSSI but 
for bookkeeping purposes by the computer they are labeled A 
through V. 
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D. INTERPRETATION OF GRAPHS 



Once all the data had been collected and the 
effectiveness for a given R and N-(- u values had been 
computed, plots of effectiveness as a function of N^ u for 
the entire range of R values were plotted for both parallel 
and counterflow heat exchangers. These plots are shown in 
Appendices H and I. 

It may be observed that for a particular value of R, the 
plot of effectiveness as a function of N tu indicates that 
the 1:5-3C exchanger outperforms the 1-4, 1-3 :2P and the 
1-5: 3P exchangers but not the 1-3 :2C exchanger. This may be 
explained by the quite apparent fact that the 1-5 :3C 
exchanger has one more counterflow pass than the 1-4 
exchanger. Indeed, this argument can be extended to the 
fact that the 1-5: 3C has one more parallel flow pass than 
the 1-3:2C. Figure 4.1. displays the graphical picture to 
support this contention. 
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Figure 4.1 Comparison of Analytical (1-4) and Computer Results 
(1-3:2C and 1-3:2P) to Computer Results (1-5:3C 
and 1-5 :3P) at R = 0.5 

46 




V. POLYNOMIAL REGRESSION 



A. DEVELOPMENT OF POLYNOMIAL EQUATIONS 

The empirical data obtained for the two 1-5 heat 
exchangers was designed to cover an extensive range of R 
values varying from 0.01 to 1.0 in increments of 0.05. As 
discussed earlier in Section IV, the data obtained at a 
specific value of R is the computer evaluated result of the 
effectiveness, for an associated N tu value. With this 
accomplished, it then becomes possible to graph separate 
curves for each of the different R values as shown in 
Appendices H and I. Through a polynomial regression 
technique, as discussed in this section, it is also possible 
to develop implicit equations for the curves with e = f(N^ u , 
R). It is also apparent from an inspection of the graphical 
representation of the empirical data in Appendices H and I, 
that the curves conform to a high degree polynomial. 

However, further analytical investigation is needed to 
ascertain the exact order of the polynomial terms. This 
investigation will not only lead to the order of the 
polynomial, but to the specific equation for each curve. 

By use of polynomial regression, the least-squares 
method can be readily extended to best fit the data to the 
m th_d e g ree f or the polynomial 

y = Aq +• A^x + A 2 X^ + ... A m x m (5.1) 
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with the error defined by 



©i — — y^ — — Aq — A^x — A2X^2 — • . . — A[jjXj_ m 

where represents the empirical data value corresponding 
to x^ , x^ being free of error. 

The objective is to minimize the sum of the squares of 
the residuals, S r , 



s = 

r 


m 2 m 2 m 2 

e f = E (D. - A Q - A l x 1 - A 2 x 2 ^ + ... A m x m ) (5.2) 


Because 


at a minimum, the partial derivatives aS r /3Ag, 


aSp/SAi 


... 3S r /aA m vanish, after taking the derivative of 


S r with 


respect to each of the coefficients of the poly- 


nomial , 


it can be seen that 


3S r 

aA o 


= 0 = -2 E (D. - A- - A,x. - A_x. 2 - ... - A x. m ) 

v l 0 li 2 i mi 


3S r 

9A 1 


- 0 = -2 E x.(D. - A_ - A, x . - A n x. 2 - ... - A x. m ) 

i v l 0 li 2i mi 


3S r 

9A 2 


= 0 = -2 E x. 2 (D. - A_ - A, x . - A 0 x. 2 - ... - A x. m ) 

i v i 0 li 2i mi 


« 


• 


• 


• 


• 


• 


as 

r 

a a 

m 


= 0 = -2 E x. m (D. - A n - A, x . - A.x. 2 - ... - A x. m ) 

l l 0 li 2 l mi 
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Then by dividing by -2 and rearranging we obtain 
AqN + Ai J + A^x-^ 2 + ... + A m E Xi m = E 

A 0 E Xi + Ai E Xi 2 + A 2 E Xi 3 + ... + A m E x-[ m+1 = E xiD-l 

Aq E Xi 2 + A^ E Xi 3 + A 2 E x-L 4 + ... + A m E x^ m+2 = E x^ 2 D^ 



A 0 E Xi 111 + Ax E x i m+1 + A 2 E x i m+2 + ... + A m E Xi 2m = E Xj^Di 

where all summations are from i=l through n. All of the 
foregoing m+1 equations are linear and have m+1 unknowns: 

Aq, Aj_, A 2 , ... A m . The coefficients of the unknowns can be 

calculated directly from the observed data. Thus, the pro- 
blem of determining a least-squares polynomial of degree m 
is equivalent to solving a system of m+1 simultaneous linear 
equations. Putting the equations in matrix form yields 



N E Xi E Xi 2 E Xi 3 . • • E xi m 




E Di 


E x-^ E x-^ 2 E X-^ 3 E x^ 4 • • • E Xi m + 1 




E 


I Xi 2 £ E E Xi 5 . . . E x 1 m+ 2 

• • * 

. • 

• • • * * 


A = 


E Xi 2 Di 


E Xi m E x i m+1 E Xi m+2 E Xi !n+3 • • • E Xi 2m 




E Xi^Di 



[Ref. 15: pp. 302-309 and Ref. 16: 468-474]. 
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From this point on, one finds that it is best to use a 
computer to assist in solving the simultaneous equations and 
this will also help alleviate any ill-conditioning that may 
otherwise occur. An existing "curvefit" program available 
through NON-IMSL [Ref. 16] was used although some modifica- 
tions were made to the original program to best accommodate 
the goals of this work. 

To determine the order of polynomial that should 
eventually be used, one increases the degree of the 
approximating polynomial as long as there is a statistically 
significant decrease in the variance o^, which is computed 
by 



2 

a 





(5.3) 



In otherwords, the selection of the optimum degree poly- 
nomial is contingent upon a decreasing variance and once the 
variance begins to increase, the degree of the polynomial 
becomes too high. For all cases, it was found that the 
e - N tu developed curves are of the 5th order. 

As shown in Figures 5.1 and 5.2 the computed values of 
effectiveness vs. N^ u for R = 0.1, 0.5 and 1.0 for both 
flow arrangements, (1-5: 3P) and (1-5: 3C), have been graphed 
and fitted by a 5th degree polynomial. Because all computed 
values for effectiveness follow a predictable trend, 
only a sample of the data covering the whole range of 
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Figure 5.1 1-5 :3C Data Fit by a 5th Order Polynomial 




Figure 5.2 1-5:3P Data Fit by a 5th Order Polynomial 
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values of R have been shown. It is clear that the graphical 
interpretation strongly backs what is known numerically 
from the polynomial regression technique. Thus, the 
relationship for e = f(N tu , R) can be found explicitly from 

e * A 5 N tu 5 + A 4 N tu 4 + A 3 N tu 3 + A 2 N tu 2 + A l N tu + A 0 (5*4) 

where the coefficients A 5 , A 4 , A 3 , A 2 , A]_, and Aq relating 
to a specific value of R are found in Tables 1 and 2 for the 
(1-5:3C) and (1-5:3P) configurations. 

In order to better utilize the coefficients found in 
Tables 1 and 2, an interactive program has been developed 
using equation (5.4) and the 1-5 heat exchanger polynomial 
approximation coefficients. This program was written for 
the IBM 3033 using the FORTRAN VS compiler (See Appendix J) . 
By entering the values for NTU, R and the type of heat 
exchanger the effectiveness values are readily obtained 
directly or by linear interpolation. 

An example of how to work with the effectiveness by 
using the Equation (5.4) follows. 

B. NUMERICAL EXAMPLE 

1. Statement of Problem 

A 1-5 :3C heat exchanger uses sea water (Cp = 4.045 
kJ/kg-°C) at an initial temperature of 20°C to cool engine 
oil, flowing at 200 kg/s from 60°C to 40°C. The sea water 
flow rate is 80 kg/s. If the overall heat transfer 
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coefficient is 40 W/m 2# C, how much surface area is required 
in the exchanger? 

2. Solution 

At an average oil temperature of 50°C, the specific 
heat is Cp = 2.004 kJ/kg-°C and an energy balance gives 

(WCpAT) sea wa ter = ) engine oil 

(80 kg/s ) (4 . 045 kJ/kg-°C)(t 2 - 20)°C 

= (200 kg/s ) (2 . 004 kJ/kg-°C) (60 - 40)°C 



t 



2 



20 + 



200(2. 004)(20) 
10(4.045) 



44. 8°C 



The number of transfer units is given by Nt u = 
UA/Cmin* To determine we must compare the products of 

mass-flow rate and specific heat for sea water and oil. 

(WCp) S w = (80 kg/s ) (4 . 045 kJ/kg-*C) 

= 323.6 kW/°C 

< wc p)oil = (200 kg/s ) (2 . 004 kJ/kg-°C) 

= 400.8 kW/°C 

Therefore C m t n = 323.6 kW/°C. 

The effectiveness is a function of R and Nt u . The 
capacity rate ratio, R, is 
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R 



323.6 

400.8 



0.8074 



C . 
min 



max 



Because C c = C m i n , the required effectiveness is 



t 2 “ t l _ 44.8 - 20 

T 1 - t 1 60 - 20 

We must now go to Table 1 for 
R = 0.8 and R = 0.85 and find 

. R = 0.8 

A 0 = 3.63419 x 10 -3 

A x = 0.93822 

A 2 = -0.67053 
A 3 = 0.27531 

A 4 = -0.60558 x 10 _1 
A 5 = 0.54393 x 10" 2 



= 0.6193 

the 1-5 :3C arrangement with 
the coefficients 

R = 0.85 



A 0 


= 


3. 82659 


X 


CO 

1 

o 

rH 


Al 


= 


0.93407 






a 2 


= 


-0.67841 






a 3 


= 


0.28109 






a 4 


= 


-0.62088 


X 


h- 4 

o 

1 

h- 4 


a 5 


= 


0.55850 


X 


10~ 2 



A trial and error solution (assuming values of N-^ u and by 
interpolating) gives a result 



N tu = 2.352 

for which 



A ” N tu Cmin/U 
A = 19.03 m 2 
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3. Observations 



The primary observation made here is that by using 
the 5th order polynomial equation with the appropriate 
coefficients found in Table 1 or 2, an accurate value for 
effectiveness can be found thus allowing one to solve for 
a wide variety of heat exchanger parameters. 

By comparing the value for effectiveness just 
computed with that of the 1-5: 3P, 1-3 :2P [Ref. 1] and 1-4 to 
1-12 even pass exchangers [Ref. 8], the 1-5: 3C has a higher 
ef f-ectiveness than all (see Table 3). 

From this observation it is evident that the 1-5 :3C 
exchanger outperforms not only the 1-4 to 1-12 even pass 
arrangements but also its counterpart, the 1-5: 3P exchanger 
for these given values by at least 0.614%. 
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TABLE 1 (cont’d) 

-5 : 3C COEFFICIENTS 
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:3P COEFFICIENTS 
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TABLE 2 (cont’d) 

-5 : 3P COEFFICIENTS 
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TABLE 3 

COMPARISON OF EXAMPLE RESULTS 
FOR Ntu = 2.352 AND R = 0.8074 



1 NUMBER OF 1 

1 TUBE PASSES 


EFFECTIVENESS 


3 : 2P 


.5990 


4 


.6155 


5 : 3C 


.6193 


5 : 3P 


. 6067 


6 


.6145 


8 


.6141 


10 


.6140 


12 


.6139 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The objective of this thesis was to provide effective- 
ness values for the 1-5 :3C and 1-5 :3P heat exchangers. 
Although these objectives were analytically unachievable, 
they became attainable through a fifth order polynomial 
approximation from numerical data. The data obtained covers 
a wide range of capacity rate ratios for values of N^ u from 
0 to 3.25. With the knowledge of the 1-5 heat exchanger 
type (1-5: 3C or 1-5: 3P), size (surface area), the overall 
heat transfer coefficient and the fluid flow rates, the 
N-£ U and capacity rate ratios may be computed and the 
exchanger effectiveness may then be determined from the 
appropriate coefficient tables or the effectiveness program 
provided . 

It is indeed obvious that the effectiveness parameter is 
most useful in performance estimation and analysis for all 
types of heat exchangers. The work done here pertaining to 
the 1-5 exchanger is very helpful in constructing a 
comprehensive picture of the performance attainable for the 
1-2 consecutively through the 1-6 for a diverse range of 
known conditions. 
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B . RECOMMENDATIONS 



With proof in this thesis, and that of O'Hare [Ref. 1], 
that an odd number of tube passes per shell has a slightly 
better effectiveness when more than half the tube fluid 
flows in counterflow to the shell fluid, it is recommended 
that a feasibility study be conducted with regard to the 
structural and thermal problems that may arise in 
manufacturing and designing such heat exchangers. 
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905 FORMAT (/ ' INPUT TOTAL HEAT TRANSFER SURFACE: 
READ (IN, 902) SURFTO 
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NTU53P COMPUTER GENERATED INPUT ANALYZER PROGRAM 
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WRITE (IOT, 904) 

904 FORMAT?/' INPUT OVERALL HEAT TRANSFER COEFFICIENT: 
READ (IN , 902 ) U 

WRITE (IOT , 905 ) 



905 FORMAT (/ ' INPUT TOTAL HEAT TRANSFER SURFACE: 
READ (IN, 902) SURFTO 

WRITE (IOT, 906) 

906 FORMAT?/' INPUT HOT SIDE INLET TEMPERATURE:' 
READ( IN , 902 ) THOTIN 
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SAMPLE OUTPUT FROM NTU53C COMPUTER INPUT ANALYZER 
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MODIFIED NTU53C PROGRAM USED TO RUN ON BATCH SYSTEM 
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Figure H . 3 



1-5: 3P Effectiveness vs, 
0.5 to 1.00 



over Range of R from 
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1-5 :3C EFFECTIVENESS VS. N GRAPHS AT VARIOUS R VALUES 
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Figure 1.2 1-5:3C Effectiveness vs. N over Range of R from 

0.01 to 0.5 
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Figure 1.3 1-5 :3C Effectiveness vs. N over Range of R from 

0.5 to 1.00 
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